Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 1 February 2008 (MN MfeK style file vl.4) 



Magnetic flares and the optical variability of the X-ray 
transient XTE J1118+480 

A. Merloni 1 , T. Di Matteo 2 * and A. C. Fabian 1 

1 Institute of Astronomy, Madingley Road, Cambridge, CBS OH A. 

2 Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138, USA. 



o 
o 
o 

(N 



(N 
> 

m 

\o 
o 
o 
o 

S3 

Oh 
6 



,H ,1 INTRODUCTION 



ABSTRACT 

The simultaneous presence of a strong quasi periodic oscillation of period ~10 seconds 
in the optical and X-ray lightcurves of the X-ray transient XTE J1118+480 suggests 
that a significant fraction of the optical flux originates from the inner part of the 
accretion flow, where most of the X-rays are produced. We present a model of magnetic 
flares in an accretion disc corona where thermal cyclo-synchrotron emission contributes 
significantly to the optical emission, while the X-rays are produced by inverse Compton 
scattering of the soft photons produced by dissipation in the underlying disc and by 
the synchrotron process itself. Given the observational constraints, we estimate the 
values for the coronal temperature, optical depth and magnetic field intensity, as well 
as the accretion rate for the source. Within our model we predict a correlation between 
optical and hard X-ray variability and an anticorrelation between optical and soft X- 
rays. We also expect optical variability on flaring timescales (~ tens of milliseconds), 
with a power density spectrum similar to the one observed in the X-ray band. Finally 
we use both the available optical/EUV/X-ray spectral energy distribution and the low 
frequency time variability to discuss limits on the inner radius of the optically thick 
disc. 

Key words: accretion, accretion discs - magnetic fields - radiation mechanisms: 
thermal - binaries: general - stars: individual: XTE Jl 18+480 
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The newly discovered transient X-ray source XTE 



J1118+480 (Remillard et al. 2000), has already shown a 
number of remarkable properties. It was discovered as a 
weak source with the RXTE All-Sky Monitor on March 29, 
2000 at high galactic latitude (b ~ 62°). Subsequent RXTE 
pointed observations revealed an energy spectrum typical of 
black hole candidates in their hard state, showing a promi- 
nent power-law component in its X-ray spectrum with a pho- 
ton index of about 1.8 up to at least 30 keV. The source was 
also observed on March 26 in hard X-rays by BATSE, up 
to 120 keV (Wilson & McCollough 2000). The 13th mag- 
nitude optical counterpart, discovered by Uemura, Kato & 
Yamaoka (2000), exhibit a spectrum fairly typical of X-ray 
Novae in outburst (e.g. Garcia et al. 2000). 

Being at such an high galactic latitude, the source suf- 
fers from low interstellar absorption, and consequently has 
also been observed in the extreme ultra violet (EUV) by the 
Extreme Ultraviolet Explorer (EUVE). The EUVE spectrum 
and fluxes derived by Hynes et al. (2000) however are ex- 
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tremely sensitive to the details of the absorbing column and 
are therefore subject to large uncertainties. 

A strong QPO feature with a frequency v ~ 0.08 Hz 
has been found in the X-ray power density spectrum (PDS) 
of the source in the first month of the outburst (March 29 
- May 4, see Revnitsev, Sunyaev & Borozdin 2000). More 
surprisingly, the optical lightcurves also show strong flicker- 
ing on timescales of a few seconds or faster and a prominent 
quasi periodic feature with a frequency in agreement with 
that of the X-ray QPO. Subsequent HST, RXTE and ASCA 
observations in the optical/UV and X-ray band respectively, 
have also shown that the QPO is drifting to higher fre- 
quencies systematically in both bands, from v ~ 0.08 Hz 
to v ~ 0.11 Hz later in the month (Haswell et al. 2000b; 
Yamaoka, Ueda & Dotani 2000; Patterson private commu- 
nication) . 

The strikingly similar variability properties observed in 
both bands suggest that at least a fraction of the optical/UV 
flux from the object should be produced in the same region 
of the accretion flow where the X-rays are produced and 
hence reflect a different aspect of the same phenomenon. 
This implies that the modulated optical/UV flux cannot be 
produced by thermal emission from the accretion disc itself; 
the optical/UV flux from a stellar mass black hole arises 
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from the outer region of the disc (radii greater than hundreds 
of Schwarzchild radii) where the X-ray emission is negligible. 

Here we show that the presence of such optical/UV vari- 
ability in the observations of XTE Jll 18+480 argues for the 
presence of significant thermal synchrotron radiation from a 
magnetically dominated corona (or possibly from an advec- 
tion dominated accretion flow, ADAF, occupying the inner 
parts of the flow). 

The only other object ever observed to show similar 
short timescale variability, in both the optical and X-ray 
bands, is the black hole candidate GX 339-4 in its hard 
state (Motch et al 1982). For the case of GX 339-4, typi- 
cally brighter than XTE J1118+480 in X-rays, it was also 
suggested that thermal synchrotron emission should provide 
a significant contribution to the optical band (Fabian et al. 
1982; Di Matteo, Celotti & Fabian 1999). 



2 CYCLO-SYNCROTRON EMISSION FROM A 
MAGNETIC CORONA 

2.1 The Magnetic field in the corona 

It is now well established that magnetic stresses in accre- 
tion discs are likely to be responsible for the transfer of 
angular momentum (see Balbus & Hawley 1998 for a re- 
view). Magneto-rotational instabilities can drive turbulence 
by amplifying the seed magnetic fields on roughly Keplerian 
timescales; the rate at which magnetic energy is built up is 
fast enough to explain the bulk of energy release in an ac- 
cretion disk as magnetic dissipation. This supports the idea 
that accretion disc coronae, the loci where the hard X-ray 
emission is produced, are highly magnetic and form by buoy- 
ancy of the strong magnetic fields amplified in the disk (e.g. 
Galeev, Rosner & Vaiana 1979). This picture is supported 
by numerical simulations (Miller & Stone 2000) which show 
that, when weak B fields are amplified via MHD turbulence 
in the disk, only a fraction ~ H/R < 1 of the energy is dis- 
sipated locally while the rest escapes and forms a strongly 
magnetized corona above the disc. 

In any such models, the magnetic field in the corona is 
likely to be strongly inhomogeneous and to dissipate energy 
in localized active regions. Following Di Matteo, Celotti & 
Fabian (1997; 1999), we assume that a significant fraction, 
/, of the accretion power, L = ml/Edd, is accumulated in 
the corona while the remaining fraction (1 — /) is dissi- 
pated internally in the disc. The magnetic field strength, 
B, in the magnetic flux tubes rising from the disk can be 
higher than the values implied by equipartition with radia- 
tion energy density. In fact, if the active region is powered 
by the release of magnetic energy, its dissipation velocity is 
~ va, the Alfven speed (see e.g. Di Matteo 1998). The en- 
ergy stored in the fie ld is therefore a fraction va/c ~ 0.1 o f 
that of the radiation (Di Matteo, Blackman & Fabian 1997). 
Based on similar arguments Haardt, Maraschi & Ghisellini 
(1994) showed that the magnetic energy is released on typ- 
ical timesclales to ~ lORb/c, where Rb is the size of the 
active region. This leads to a magnetic field strength, 



sion in units of R B = 2GM/c 2 , the Schwarzschild radius for 
a given central mass M — m Mq . This implies a value 
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for the estimated magnetic field in the inner part of the 
accretion flow. 



2.2 Synchrotron emission 

In the hot (T ~ 50 — 200 keV), highly magnetized, coro- 
nal plasma the majority of the electrons are thought to be 
thermal and are expected to produce synchrotron emission, 
and Compton scatter the softer photons produced both by 
dissipation in the underlying accretion disk and by the syn- 
chrotron processes themselves. 

In a thermal plasma, optically-thin synchrotron emis- 
sion rises steeply with decreasing frequency. Under most cir- 
cumstances the emission becomes self-absorbed and gives 
rise to a black-body spectrum below a critical frequency: 
v, a 2.7 x 10 m t - 08 ( t ^ e )°- 96 ( 1 ^ 5 ) - 91 Hz (see Wardziriski 
& Zdziarski 2000), where T is the coronal temperature and 
r the scattering optical depth. Above this frequency, which 
is right in the optical/UV band for typical parameter of X- 
ray transient sources and magnetic fields B ~ 10 6 — 10 7 G, 
the cyclo-synchrotron emission decays exponentially as ex- 
pected from a thermal plasma, due to the superposition of 
cyclotron harmonics. 

In an active region above an accretion disc, the rel- 
ative importance of the cyclo-synchrotron emission (and 
its Comptonization) as opposed to the Comptonization of 
the external photon field has to be estimated by compar- 
ing the relevant timescales for the processes. Such com- 
parisons have shown that, as the synchrotron emission is 
heavily self-absorbed, the internal energy density of the 
synchrotron radiation C/ ra d,int < B 2 /8tv < Crad.ext, where 
£/rad,cxt = (1 — /) L / A-KcR\ isc is the 'external' energy den- 
sity of the radiation emerging from the disc. Consequently 
(and in accordance to standard models) most of the X-ray 
emission in the soft and hard states of GBH is usually due to 
Comptonization of the external photon fields (and not of the 
synchrotron component; Di Matteo, Celotti & Fabian 1997; 
1999; Wardziriski & Zdziarski 2000). Nevertheless, in the 
hard states of GBH the disk blackbody component, which 
dominates the X-ray emission in the soft state, often dis- 
appears or cools down significantly, and the synchrotron 
component is likely to contribute at least in the optical/UV 
band. 

Here we show that thermal synchrotron emission can 
be an important radiation process in the case of XTE 
J1118+480. 



3 APPLICATION TO XTE J1118+480 
3.1 Spectral modeling 

Hubble Space Telescope observations of XTE J1118+480 
(Haswell et al. 2000a) were made on 2000 April 8, at the 
end of the rising phase of the outburst. The data indicate 
a flux density at 1500A of 3 x 10" 13 ergs cm~ 2 s" 1 A~ 1 and 



where N is the number of active regions and rb their dimen- a spectrum between 1200 and 8000 A slightly flatter than 
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On the other hand, unfiltered CCD fast pho- 



tometry ( |Haswell et al. 20001: ; Patterson 2000 ) has revealed 
optical flickering of ±0.2 mag on timescales of seconds, which 
correspond to ~ 20 per cent of the optical flux. As the opti- 
cal variability is strongly modulated by a QPO with a fre- 
quency comparable (if not the same) to that of the X-ray 
one, it has to be produced where most of the X-rays are 
produced, namely in the inner part of the accretion flow. 

We assume that this rapid variability is due to self ab- 
sorbed cyclo-synchrotron (CS) emission. To explicitly show 
the various parameter dependences of the synchrotr on emis- 



sion we write the flux at th e peak frequency v c as (Di Mat- 
teo, Celotti & Fabian 19971) 
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where d is the distance of the source in kiloparsecs, 9 — 
kT/m e c 2 the dimensionless temperature, U Ta d = U ra d,int + 
C^rad.ext and the ratio of the relevant timescales for the 
cyclo-synchrotron and inverse Compton emission determines 
which process dominates in an active (see Section 2.2). Un- 
der the assumption that at least one fourth of the flux 
at v = v c = 10 15 Hz is due to CS, we deduce that the 
intrinsic dissipation in the accretion disc has to be fairly 
small (m(l — /) ~ 10 -3 ) implying a magnetic field intensity 
B ~ 2 x 10 B G. 

We model the spectrum of the source following the work 
of Di Matteo, Celotti & Fabian (1999). The main feature of 
the model is that it considers reprocessing of coronal ra- 
diation in the accretion disc not only according to the ac- 
tive blob size but also according to their height above the 
disc. The temperature and the optical depth in the corona 
are constrained from the observed slope of the power-law in 
the X-ray spectrum (a ~ 0.8 oc (r#) -1 ; e.g. Wardzihski & 
Zdziarski 2000) . The spectrum is calculated self-consistently 
for every annulus of radius R = rRs and width dR <C R and 
integrated over radius from r = 3 to r = 1000. We take into 
account all the relevant radiative processes and rescale them 
opportunely according to the ratios of their typical cooling 
timescales which are obtained by integrating each spectral 
component (more accurately than in Eq. 3) 

In all cases we postulate that the accretion disc extends 
down to the innermost stable orbit of a non-rotating black 
hole (Rin = 3Rs), so that both the optical and X-rays are 
produced where most of the energy is dissipated. The num- 
ber of active regions at any given annulus is calculated as in 
Haardt, Maraschi & Ghisellini (1994), following their discus- 
sion on the timescales over which the magnetic field is am- 



plified in the disc (t z 



ixcp) and released in the corona 



(i re i ~ to, see section 2.1). We obtain dN = 9.5r~ 3 ^ 2 dr, so 
that the total number of flares active at any time N ~ 10, 
as usually assumed from variability arguments. The mass is 
fixed to 10 solar masses. 

In Figure 1 we show an illustration of two model spec- 
tra for the recent optical and X-ray observations of XTE 
Jl 118+480. They are not to be regarded as detailed fits to 



the data, rather as indicative of the physical properties of 
the source implied by the model. In particular, we do not at- 
tempt to fit the optical continuum with a realistic model for 
the outermost parts of the accretion disc, our main interest 
being to put constraints on the nature of the inner accretion 
flow. The relatively high optical-to-X-ray flux ratio for the 
source imposes the major constraints on the model. In order 
for both the synchrotron component to be important in the 
optical band and the Comptonized flux not to exceed the 
observed limit in the hard X-ray band (7 ra d, ext, the external 
radiation energy density intercepting the active region (see 
Eq. 3) has to be lower than (1 — f)L/(iirR 2 lia( ,c). 

The first spectrum (case a, solid lines) is obtained as- 
suming static dissipation regions at fixed height. For n, = 3 
we have to allow the reconnection sites to be above the ac- 
cretion disc at a height of ~ 4 — 5r^. The need for such 
aspect ratio of the active region is just an indication that 
the soft radiation field seen by the flaring region has been 
reduced. 

As expected from the simple scaling arguments of 
Eqs. (2) and (3), we obtain a relatively low accretion rate 
(m ~ 0.01) and a high value of the fraction of the power 
dissipated in the corona (/ ~ 0.97), both consistent with 
the source being in its hard state. However, in this geom- 
etry / is not a relevant parameter because the hard X-ray 
radiation illuminates the disk and a substantial fraction of 
it (^ 0.5, depending on the disc albedo) is reprocessed and 
thermalized giving rise to a hotter blackbody-like compo- 
nent roughly similar to that obtained for lower values of /. 

In Figure 1 (case b, dashed lines) we illustrate an al- 
ternative way to reduce the contribution from this thermal 
component in the soft X-ray band and enhance the the cyclo- 
synchrotron emission in the optical/UV band. Following Be- 
loborodov (1999), we allow for bulk relativistic motion of 
the emitting coronal plasma away from the disc. Due to 
Doppler boosting the radiation from a reconnection site to- 
wards the disc is reduced. This is equivalent, for our purpose, 
to the assumption that the coronal emission is anisotropic: 
for v/c ~ 0.3, roughly 90 per cent of it is emitted upwards 
and only the remaining ~ 10 per cent impinges on the un- 
derlying disc. 

Finally, the values of temperature and optical depth of 
the active regions we derive are, respectively, ~ 0.3 and 
t ~ 0.6 and are almost equal in the two cases. The distance 
we obtain is of the order of D = 0.4 kpc, as expected for a 
high latitude galactic source. 

3.2 The inner disc 

Hynes ey al. (2000) have recently reported on EUV observa- 
tions of XTE J1118+480 and have shown that the relatively 
low EUV flux may be inconsistent with an optically thick 
disc extending down to radii r[ n < 1000. However EUV ob- 
servations are extremely sensitive to the assumed absorption 
and the derived limits on r- m highly uncertain. This is further 
emphasized by recent ASCA observations which appear to 
show a slight soft excess below 2 keV (Yamaoka et al. 2000) 
well fitted by a multicolour disc of temperature 0.2 ± 0.1 
keV. Because of the inconsistency between the ASCA and 
EUVE measurements we did not use the EUVE data as a 
further constraint for our model. 

However, in accordance with the requirements discussed 
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.2.1 The 0.1 Hz QPO 
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Figure 1. Representative model predictions for the spectral en- 
ergy distribution of XTE J1118+480. The points represent op- 
tical (V band) (Uemura et al. 2000), near UV (Haswell ct al. 
2000a) and X-ray (Remillard et al. 2000, Yamaoka ct al. 2000) 
observations of the source during its outburst. The two thicker 
lines represent models with (case a, solid line) and without (case 
b, dotted lines) bulk relativistic motion of the active regions. 
The parameters are very similar for the two cases: for the ac- 
tive region parameters we obtain, respectively r = 0.61, 8 = 0.3, 
rj, = 2.9, h = 5.3rb (case a) and r = 0.56, 8 = 0.33, rj, = 3.2, 
h = 4rt, (case b). We have assumed a mass of 10 Mq, and we 
estimate a distance of 0.4 kpc; the accretion flow parameters are: 
m = 0.01 and / = 0.97. The self absorbed cyclo-synchrotron 
emission and its Comptonization arc plotted (thinner lines). The 
optically thick part significantly contributes to the optical emis- 
sion, together with the expected (non- varying) emission from the 
outermost parts of the accretion disc, while the X-ray are mainly 
produced by Comptonization of disc photons. 



above, most objects observed in their hard/low state often 
show either no or very little evidence for soft blackbody 
emission or strong reflection features (i.e., the backscattered 
emission from the putative accretion disk). Based on this, 
it has therefore been argued that, in the hard state, geo- 
metrically thin discs may not extend down to the inner- 
most stable orbit but are truncated at tens or hundreds 
of Schwarzchild radii. (Gierlihsky et al. 1997; Poutanen & 
Coppi 1998; Zdziarski et al. 1998; Done & Zycki 1998; Esin 
et al. 1997; Wilms et al. 1999 and references therein). 

If the disc does not extends down to the innermost sta- 
ble orbit, then [7 ra d, C xt is suppressed and the relative impor- 
tance of the synchrotron component (see Eq. 3) increases, as 
required by the data. Therefore, a central advection domi- 
nated (ADAF) component of the flow extending out to hun- 
dreds or thousands of Schwarzschild radii and surrounded by 
an optically thick accretion disc may also fit the spectrum 
reasonably well (see e.g. Esin et al. 1998). The information 
provided by temporal studies are then crucial for constrain- 
ing the geometry of the inner disk, in particular optical- 
X-rays cross-correlation and time lag analysis for the QPO 
component. 



Psaltis, Belloni & van der Klis (1999) have suggested that 
the frequencies of the periodic features in the PDS of black 
hole candidates (which can appear more or less broadened) 
follow a quite tight correlation that encompasses also neu- 
tron star sources, and therefore has to be related to the 
properties of the accretion flow. Recently Psaltis & Norman 
(2000) have proposed a model in which the QPOs are pro- 
duced in a narrow ring in the accretion disc where a discon- 
tinuity in some properties of the flow occurs. The frequen- 
cies themselves are essentially determined by the relativistic 
proper frequencies of the system (Stella, Vietri & Morsink 
1999; Merloni et al. 1999). Assuming such a model, the data 
can be used to constrain the position of the flow discontinu- 
ity (Di Matteo & Psaltis 1999). 

The published PDS of XTE J1118+480 strongly resem- 
bles that of GX 339-4 in the hard state (see Nowak 2000 for 
a discussion of the various variability components). Apart 
from the strongly peaked QPO at 0.08 Hz, it exhibits a 
broader feature around 1 Hz. Plotted one against the other, 
these frequencies fall remarkably well on the correlation pre- 
sented by Psaltis, Belloni & van der Klis (1999). This would 
suggest that the low frequency QPO is analogous to the 
HBO in neutron stars and is related to twice the nodal pre- 
cession general-relativistic frequency of a perturbed orbit 
around the compact object. The value of the transition ra- 
dius, where the modulation is produced, depends on the spin 
of the black hole, but is limited in the range 25 ^ Rt /Rs 6 
for a 10 Mq black hole with dimensionless angular momen- 
tum 0.99 a 0.01. This provides an upper bound for the 
inner extent of the geometrically thin accretion disc, which is 
in conflict with the limits derived from EUVE observations 
(Hynes et al. 2000) and significantly reduces the relevance 
of a central ADAF component, favoring instead models dis- 
cussed in Section 3.1. We note, however, that if the observed 
QPO is associated with the HBO frequency it is not clear 
why the expected modulation at the keplerian frequency at 
the same radius (with frequency i* cp ~ 30 Hz) is not di- 
rectly observed, even though it should have a higher quality 
factor (although, see e.g. Nowak 2000). 

Alternatively, we can ignore the very broad feature in 
the PDS at ~ 1 Hz and consider the sharp low frequency 
QPO as due to Keplerian modulation of the flow far away 
from the source. This would imply a transition radius R t ~ 
500-Rs- Although this interpretation for the observed vari- 
ability is less plausible, it may be in better agreement with 
the implied EUV fluxes and possibly with the presence of a 
central ADAF component. However, in this case it remains 
to be explained how the modulation propagates inward in 
the advective flow without loss of coherence, given that, in 
the ADAF model, the bulk of the optical and X-ray are pro- 
duced in the inner part of the flow, where the magnetic field 

1/2 . 1/2 -5/4 



is high enough (Badaf ~ a o\ /2 m u)' ^ ™osn r z°' ^ G) to 
power the cyclo-synchrotron emission, which is then Com- 
ponized to give rise to the high energy spectrum (see e. g. 
Narayan & Yi 1995; Quataert & Narayan 1999). 



4 DISCUSSION 

We have applied a model for the emission from a highly 
magnetic, structured corona to the optical and X-ray ob- 
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servations of the newly discovered X-ray transient XTE 
J1118+480. 

This source is unique for a number of reasons. It is lo- 
cated at high galactic latitude and has a very high optical- 
to- X-ray flux ratio. The source is observed in the typical 
black hole candidates hard state (photon index V ~ 1.8). 
Given the optical/UV to X-ray flux ratio we have derived 
constraints for the approximate size, optical depth and mag- 
netic field strengths of the coronal active regions. 

The main point of our work is that the simultaneous 
presence of a strong quasi periodic feature in the optical 
and X-ray lightcurves clearly suggests that the fluxes in the 
two bands both originate from the same region in the inner 
part of the accretion flow. Self-absorbed cyclo-synchrotron 
emission is the natural candidate to explain the optical vari- 
ability. Such emission is expected in any magnetic corona 
model, and the inferred magnetic field value (B ~ 2 x 10 6 
G) is the one predicted to arise when the source is in the 
hard state (low value of m and high value of /). 

Our model is not unique. The relatively high optical- 
to-X-ray flux ratio observed can be explained either by a 
strong, structured magnetic corona with a relatvely high 
scaleheight, or by a central ADAF surrounded by an op- 
tically thick accretion disc. We have shown that even if 
cyclo-synchrotron emission plays an important role in the 
optical band in both cases, in our model the X-ray photons 
are mainly produced via inverse Compton scattering of the 
soft disc photons, while in the ADAF model the CS emis- 
sion itself acts as source field. We gave an example of how 
variability data can be used to further discriminate between 
the models. 

Independently of the dominant source of soft photons 
for Comptonization, we expect that any (small) variation 
in the physical parameters characterizing the active coronal 
regions (either in temperature or optical depths) driven by 
a flow discontinuity in the disk should modulate the emis- 
sion simultaneously in the optical (via the CS emission) and 
in the X-ray band (via inverse Compton emission). In par- 
ticular, any increase (at a level of 20 per cent or so, as re- 
quired) in temperature or density in the corona will cause 
an increase in the cyclo-synchrotron component at optical 
frequencies, while at the same time it will make the X-ray 
spectrum harder. From any such model we then expect a 
correla tion heiv 



jetweeri the optical and haid X-ra^y variability 
but an hntioorrolation between optical and ooft X rayo (thia 
has been observed for the case of GX 339-4; see Motch et al. 
1982). Our model also predicts optical variability on flaring 
timescales (~ tens of milliseconds), with a PDS similar to 
the one observed in X-rays. Furthermore the variable optical 
component is expected to drop out quickly at low frequencies 
v < v c due to the fast decline of the synchrotron emission in 
the Rayleigh- Jeans regime. Time-resolved spectroscopy ob- 
servations of such behaviour would support this model and 
allow us to place strong constraints on the magnetic field 
strength in the corona. 
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